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ABSTRACT 

Kilonovae represent an important electromagnetic counterpart for compact binary mergers, 
which could become the most commonly detected gravitational wave (GW) source. Follow¬ 
up observations, triggered by GW events, of kilonovae are nevertheless difficult due to poor 
localization by GW detectors and due to their faint near-infrared peak emission that has limited 
observational capability. We show that the Near-Infrared Camera (NIRCam) on the James Webb 
Space Telescope (JWST) will be able to detect kilonovae within the relevant GW-detection range 
of ~ 200 Mpc in short (,$ 12-second) exposure times for a week following the merger. Despite this 
sensitivity, a kilonova search fully covering a fiducial localized area of 10 deg 2 will not be viable 
with NIRCam due to its limited field of view. However, targeted surveys may be developed to 
optimize the likelihood of discovering kilonovae efficiently within limited observing time. We 
estimate that a survey of 10 deg 2 focused on galaxies within 200 Mpc would require about 13 
hours, dominated by overhead times; a survey further focused on galaxies exhibiting high star- 
formation rates would require ~5 hours. The characteristic time may be reduced to as little as 
~4 hours, without compromising the likelihood of detecting kilonovae, by surveying sky areas 
associated with 50%, rather than 90%, confidence regions of 3 GW events, rather than a single 
event. On detection and identification of a kilonova, a limited number of NIRCam follow-up 
observations could constrain the properties of matter ejected by the binary and the equation of 
state of dense nuclear matter. 

Subject headings: gravitational waves — infrared: general — methods: observational 

With the upcoming completion of advanced 
GW detectors, compact binary mergers are ex¬ 
pected to be observed via GWs; such mergers 
will likely be the first sources from which GWs 
are directly detected. The construction of Ad¬ 
vanced LIGO (Harry & LIGO Scientific Collab¬ 
oration 2010) will finish in 2015, with the de¬ 
tectors gradually reaching their design sensitiv¬ 
ity by 2019 (LIGO Scientific Collaboration et al. 
2013). Advanced Virgo (The Virgo Collabora¬ 
tion 2009) has a similar schedule with about a 
year delay. On reaching its design sensitivity, 
the Advanced LIGO-Virgo network could be able 
to detect a binary neutron star merger out to 
~ 450 Mpc under favorable conditions, or, on av- 


1. Introduction 

Compact binary mergers represent one of the 
most actively studied astrophysical phenomena. 
The mergers of neutron stars and stellar-mass 
black holes are promising targets for gravitational- 
wave (GW) observations (Abadie et al. 2010). 
They can also drive relativistic outflows, giving 
rise to some of the liighest-energy emission of 
gamma rays, neutrinos and possibly cosmic rays. 
The study of binary mergers can help us to further 
understand high-energy physics, the properties of 
matter at nuclear densities, stellar evolution, the 
abundance of heavy elements in the universe, as 
well as galaxy formation. 
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erage, out to ~ 200 Mpc (Bartos et al. 2013; LIGO 
Scientific Collaboration et al. 2013). The ex¬ 
pected detection rate of binary neutron star merg¬ 
ers with Advanced LIGO-Virgo’s design sensitiv¬ 
ity is 0.2 — 200 per year (LIGO Scientific Collab¬ 
oration et al. 2013). Additionally, neutron star- 
black hole mergers will be detectable out to about 
a factor of two farther, corresponding to compa¬ 
rable detection rates (Abadie et al. 2010). The 
construction of additional GW observatories, such 
as KAGRA (Somiya 2012) and LIGO-India (Iyer 
et al. 2011), will further increase the number of 
detectable sources. 

To maximize the sensitivity and scientific re¬ 
turn of GW observation campaigns, GW candi¬ 
dates will be followed up by other instruments 
to find electromagnetic or neutrino counterparts. 
There are significant theoretical (e.g., Metzger & 
Berger 2012; Nakar & Piran 2011; Abadie et al. 
2010) and observational (e.g., Kanner et al. 2008; 
Abadie et al. 2012; Smith et al. 2013; Bartos et al. 
2014; Aasi et al. 2014) efforts to identify promising 
counterparts. 

Kilonovae are produced during the merger of 
a binary neutron star or a neutron star-black 
hole system (Li & Paczynski 1998; Kulkarni 2005; 
Rosswog 2005; Kulkarni & Kasliwal 2009; Met¬ 
zger et al. 2010). If the neutron star is dis¬ 
rupted during the merger, the ejected neutron- 
rich matter produces heavy r-process elements, the 
radioactive decay of which powers the emission. 
Due to the large opacity at optical wavelengths 
of the formed r-process elements, kilonovae peak 
in the near-infrared with expected luminosities of 
~ 10 40 — 10 41 ergs _1 , and last for over a week 
Barnes & Kasen 2013. Nevertheless, this standard 
theoretical estimate has important uncertainties 
that may affect detectability. There is a signifi¬ 
cant ongoing theoretical effort to obtain a more 
complete understanding of the emission process, 
which may inform observations in the future (e.g., 
Kasen et al. 2015a; Fernandez et al. 2015). 

Kilonovae are promising electromagnetic coun¬ 
terparts of binary mergers because (i) the emis¬ 
sion is isotropic; therefore, the number of observ¬ 
able mergers is not limited by beaming; (ii) the 
week-long emission period allows sufficient time 
for follow-up observations; and (iii) once identi¬ 
fied, source location can be accurately recovered, 
allowing for the identification of the host environ¬ 


ment and the search for counterparts in other elec¬ 
tromagnetic regimes. For comparison, some other 
electromagnetic counterparts, such as gamma-ray 
bursts and X-ray afterglows, are highly beamed, 
reducing the number of observable mergers. Op¬ 
tical and radio afterglow is also emitted off-axis; 
however, detection at large viewing angles is dif¬ 
ficult (van Eerten & MacFadyen 2011). Gamma- 
ray and X-ray emission are also short duration, 
which can also be a significant limitation as it al¬ 
lows less time for follow-up observations. On the 
other hand, the observability of kilonovae is cur¬ 
rently limited by the lack of sufficiently sensitive 
survey instruments in the near-infrared band that 
can provide coverage over tens of square degrees, 
the typical area within which GW events will be 
localized by the Advanced LIGO-Virgo network 
(LIGO Scientific Collaboration et al. 2013). The 
sole kilonova observation so far took advantage of 
the precisely reconstructed source direction (Tan- 
vir et al. 2013; Berger et al. 2013). 

The James Webb Space Telescope (JWST) is a 
planned, highly sensitive infrared space telescope 
with an expected launch in 2018. In this paper, 
we show that its Near Infrared Camera (NIRCam; 
Horner & Rieke 2004), with a spectral range of 0.6- 
5 fin i, is well suited for quickly detecting kilono¬ 
vae following GW triggers from Advanced LIGO- 
Virgo. 

In Section 2, we determine the sensitivity of 
NIRCam to detecting kilonova emission within its 
field of view. In Section 3, we motivate the need 
for observation strategies that target nearby galax¬ 
ies to host kilonovae; in Section 4, we discuss ex¬ 
amples of such targeted surveys to identify kilo¬ 
novae. The possible role of other astrophysical 
sources being misidentified as kilonovae are pre¬ 
sented in Section 5. In Section 6, we investigate 
how NIRCam could probe the kilonova emission 
model and source parameters. We summarize our 
findings in Section 7. 

2. NIRCam sensitivity to kilonovae 

We first calculate the JWST/NIRCam sensitiv¬ 
ity in detecting kilonovae. We determine which 
NIRCam filter is the most sensitive for the ex¬ 
pected emission spectrum, and derive the integra¬ 
tion time necessary for a 10 a detection. Only then 
may we develop viable observing strategies with 
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the greatest likelihood of enabling kilonovae iden¬ 
tification. 

We adopt the kilonova emission model obtained 
by Barnes & Kasen (2013) using time depen¬ 
dent, multi-wavelength radiative transport calcu¬ 
lations. Barnes & Kasen (2013) simulated a range 
of emission parameters, in particular ejecta masses 
~ 10~ 3 — 10 _1 Mg and characteristic ejection ve¬ 
locities (3 = v/c ~ 0.1 — 0.3. These ranges seem 
to cover the expected and observed kilonova emis¬ 
sion parameters (Tanvir et al. 2013; Berger et al. 
2013). 

To determine the observed kilonova flux at 
Earth, we assume a source luminosity distance 
of 200 Mpc, which represents the average reach 
of the Advanced LIGO-Virgo network at design 
sensitivity to binary neutron star mergers. While 
some sources may occur even closer or as far as 
~ 450 Mpc under the most favorable conditions, 
our choice of 200 Mpc provides a good estimate of 
the potential limitations of a follow-up search with 
NIRCam. 

With kilonova spectra expected to peak in the 
near-infrared, we calculated the minimum integra¬ 
tion time sufficient to detect a kilonova for each 
of NIRCam’s wide-band filters. Preliminary filter 
sensitivity curves were taken from JWST’s web¬ 
page 1 . The two wide-band filters with the longest 
wavelengths [F356W, F444W] were excluded from 
the analysis because the simulations of Barnes 
& Kasen (2013) did not fully cover these bands, 
leaving us with filters F070W, F090W, F115W, 
F150W, F200W and F277W 2 . Nevertheless, as we 
will see, F200W and F277W are already suffi¬ 
ciently sensitive, and F277W is less sensitive than 
F200W soon after the binary merger, indicating 
that F356W and F444W are likely not more opti¬ 
mal for kilonova detection. 

Emulating the JWST NIRCam Exposure Time 
Calculator 3 , we determined the minimum time 
needed for the different filters to obtain a 10 a de¬ 
tection as a function of the time after a binary 
merger at 200 Mpc for the range of kilonova emis- 


http://www.stsci.edu/jwst/instruments/nircam/ 
instrumentdesign/filters 

For each filter, the number in the filter name represents 100 
times the central wavelength in microns. For example, the 
central wavelength of F356W is ~3.56 p,m.) 

'http://jwstetc.stsci.edu/etc/input/nircam/imaging/ 



Fig. 1.— Integration time for JWST NIRCam 
necessary to detect a kilonova (10 -2 M 0 ejected 
mass, (3 = 0.2) with 10 cr significance as a function 
of time after a binary merger at 200 Mpc, for dif¬ 
ferent NIRCam filters. The results indicate that 
the F200W and F277W filters are well suited to 
observe kilonovae, depending on how quickly the 
observations are obtained following the merger. 

sion parameters. In particular, for each filter, we 
first calculated the signal-to-noise (SNR) of a kilo¬ 
nova as a function of the (i) total exposure time 
and (ii) time-varying source flux, determined by 
integrating the kilonova flux density weighed with 
the filter transmittance. We then selected that ex¬ 
posure time associated with 10er, and repeated the 
process for a range of elapsed time since merger. 

Results for a kilonova with ICC 2 M 0 ejected 
mass and (3 = 0.2 characteristic velocity are shown 
in Fig. 1 as an illustration of our findings, in gen¬ 
eral. One can see that integration times of only (i) 
5-10 seconds are necessary in F277W if the kilo¬ 
nova is seen within 5 days; and, (ii) 10-12 seconds 
are necessary in F200W if it is caught within a 
week. JWST/NIRCam will have unprecedented 
sensitivity to detect kilonovae. 

3. Full survey of GW trigger events 

While NIRCam has unprecedented sensitivity, 
there are several limiting factors, relevant for a 
kilonova search aiming to cover a large, contigu¬ 
ous region of sky: (i) NIRCam has a minimum 
integration time of 10.6 seconds due to the read¬ 
out time of the detectors in full-array mode; and 
(ii) the 5" gaps between the four detectors of each 
module and the 50" gap between the two modules, 
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which necessitate positional shifts of less than half 
the 2.2' x 4.4' held of view for full coverage in the 
short wavelength channel in survey mode. These 
positional shifts would also serve to provide, for 
each point imaged, 2-3 dithered exposures, which 
may be used to remove image artifacts, such as 
bad pixels and cosmic rays. To uniformly cover 
a fiducial localized sky area of 10 deg 2 associated 
with a GW source, at a confidence level of 90%, 
these two factors would necessitate ~ 50 hours of 
total exposure time, not including overheads , mak¬ 
ing it unfeasible to cover the full region associated 
with a GW trigger event. 

4. Targeted surveys 

Several approaches may be used to develop a vi¬ 
able NIRCarn survey, triggered by a GW source, to 
search for a kilonova most efficiently and with the 
greatest chance of success. The survey should be 
triggered only by the most promising GW sources: 
those with high signal-to-noise ratios and well lo¬ 
calized areas. While we consider a fiducial sky 
area of 10 deg 2 , a substantial fraction of GW trig¬ 
ger events will have larger reconstructed sky areas 
(e.g., LIGO Scientific Collaboration et al. 2013; 
Singer et al. 2014) and therefore require propor¬ 
tionally longer NIRCarn observing times. In this 
section, we discuss approaches that yield targeted 
surveys in order to reduce the required NIRCarn 
observing time. 

Instead of covering the full 10 deg 2 area, one 
approach is to focus the follow-up search toward 
known galaxies (Kopparapu et al. 2008; Kanner 
et al. 2008; Nissanke et al. 2013; Hanna et al. 
2014; Bartos et al. 2015) within ~ 200 Mpc, taking 
advantage of the expectation that binary merg¬ 
ers occur in or near galaxies (Fong et al. 2010). 
Following Nissanke et al. (2013), the number den¬ 
sity of galaxies within 200 Mpc is estimated to be 
~8deg -2 ; therefore, one would need to follow up 
~ 80 galaxies for a fiducial GW sky area of 10 deg 2 . 

Since most galaxies can be covered well within 
a 400 x 400 pixel 2 (13" x 13") region (e.g., Fong 
et al. 2010), NIRCarn can be used in a subarray 
imaging mode (e.g., Beichman et al. 2014) for such 
a survey of galaxies. In this mode, exposure times 
as short as ~2 seconds are possible. In principle, 
this minimum exposure could be used with F277W 
for each of 3 dithered images, to mitigate image 


artifacts, of a galaxy if observed within 3 days of 
the merger. For our purposes, we consider instead 
3 dithered 4-second exposures, enabling the use of 
F200W for up to a week after the merger. For 
the 80 galaxies, assuming the most time-intensive 
case of only one galaxy present in a field, the total 
exposure time would be 16 minutes, not including 
overheads. While such an approach represents a 
significant improvement over full coverage of the 
GW event, the time associated with overheads is 
expected to be significant. 

The overheads for such a NIRCarn imaging sur¬ 
vey may be attributed to a number of sources, 
summarized in Table 1 (Gordon et al. 2012). The 
most significant contributors to overhead times 
would be the time for JWST to slew and acquire 
guide stars. The standard assumption for slew 
time to a new field is 30 minutes, which we adopt 
for the slew to the first galaxy following a GW 
event. Slew times to subsequent galaxies associ¬ 
ated with that event depend on the angular sep¬ 
aration between the galaxies. While JWST slews 
at a nominal rate of at least 90° per hour, the 
slew rates are non-linear with distance due to in¬ 
ertia. Shorter slew distances are associated with 
slower effective rates. A survey of 80 galaxies over 
a 10 deg 2 region suggests an average separation of 
~ 21' between nearby galaxies. According to ex¬ 
pectations for slew rates as a function of slew dis¬ 
tance (Gardner et al. 2010), our typical slew time 
would be no more than ~ 3.6 minutes 4 to point 
from one galaxy to the next. This typical slew 
is sufficiently large that each field will require a 
different guide star and 4 minutes for its acquisi¬ 
tion. Accounting for the move to the first galaxy, 
moves to subsequent galaxies, and guide star ac¬ 
quisitions for all fields, the total overhead time as¬ 
sociated with slewing and guide stars is 10.6 hours. 
The next most significant contributor to overhead 
times would be the time associated with the small 
(e.g., < 2") dithering, which totals 1.3 hours for all 


4 We discuss this slew time as an upper limit since it repre¬ 
sents an extrapolation from 2.5 minutes for a 5' slew and 3.3 
minutes for a 17 ' slew (Gardner et al. 2010). Since greater 
slews are more efficient, the actual slew time may be less 
than that extrapolated. Note that the 2.5 minutes for the 
5 ' slew with the 4 minutes for the guide star acquisition is 
comparable to mission requirement MR-180 (Bogenberger 
2007), reflected also in the operations requirement MO-442 
(Jordan 2014), that a 4.7' slew, including guide star acqui¬ 
sition, be accomplished in 8 minutes or less. 
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Table 1: Targeted Surveys of GW Trigger Event 



90% confidence region 

50% confidence region 


All 

Ha 

All 

Ha 


galaxies 

galaxies 

galaxies 

galaxies 

Fiducial sky area [cleg 2 ] 

10 

10 

2 

2 

Number of Galaxies 

80 

30 

16 

6 

Slew Time to First Galaxy [min] 

30.0 

30.0 

30.0 

30.0 

Slew Time to Each Subsequent Galaxy [min] 

3.6 

4.5 

3.6 

4.5 

Guide Star Acquisition Time for Each Galaxy [min] 

4.0 

4.0 

4.0 

4.0 

Dithering Time for Each Galaxy 0 [min] 

1.0 

1.0 

1.0 

1.0 

Other Overhead Time 6 [min] 

25.1 

10.1 

5.9 

2.9 


Total Overhead Time [hr] 

12.3 

5.3 

8.5 C 

4.3 C 

Total Exposure Time [min] 

16.0 

6.0 

9.6 C 

3.6 C 

Total Time [hr] 

12.6 

5.4 

00 

o 

4.4 C 


a Two dither motions, each requiring 30 seconds, in order to obtain three dithered images for each galaxy. 

b Includes 65 seconds for filter move and detector configuration, necessary only once, and an adopted detector deadtime of 
18 seconds for each galaxy. The detector deadtime for our subarray observations was assumed to be 6 seconds, based on 
discussion in Gordon et al. (2012), for each of the three images of a galaxy. 

c Total times listed for targeted survey of 50% confidence regions are those associated with three GW trigger events for direct 
comparison to total times for targeted survey of the 90% confidence region of a single event, yielding a similar likelihood of 
kilonova detection. 


fields. Finally, other overheads include moving the 
filter wheel, detector configuration and deadtime, 
which total only 0.4 hour. With the full overhead 
time of 12.3 hours, the total time for this targeted 
survey of 80 galaxies in 10 deg 2 , associated with a 
GW trigger event, is 12.6 hours. 

A second approach is to further focus the survey 
to target only those nearby galaxies most likely 
to host kilonovae. For example, the star forma¬ 
tion rate is correlated with the rate of compact 
binary mergers (Leibler & Berger 2010); Ha: is an 
indicator of the star formation rate and therefore 
may be used to optimize the probability of dis¬ 
covering kilonovae efficiently within a limited ob¬ 
serving time. Specifically, Ha mapping may be 
used to identify those galaxies responsible for 90% 
of star formation and for 50% of the mass (Met¬ 
zger et al. 2013; Bartos et al. 2015). Based on the 
galaxy stellar mass function for the nearest galax¬ 
ies (Baldry et al. 2012; Moustakas et al. 2013), 
with the component due to star-forming galaxies 
normalized such that they account for 60% of mass 
of the nearest galaxies (Metzger et al. 2013), these 
Ha-detected galaxies would represent about 40% 


of the total number of galaxies. Thus, with only 
~ 30 galaxies most likely to host kilonovae within 
200Mpc over 10deg 2 , and typical separations of 
~35', the total overhead associated with slew¬ 
ing and guide star acquisition is further decreased 
to no more than 4.7 hours. With the additional 
overheads, the total time for this targeted survey 
would be 5.4 hours. 

We note that a focus on Ha serves to solve an¬ 
other problem: a comprehensive galaxy catalog 
complete to 200 Mpc may not be ready by 2019, 
when Advanced LIGO detectors will reach their 
design sensitivities. For reference, recent galaxy 
catalogs are estimated to be only ~ 60% complete 
to 100Mpc (e.g., White et al. 2011). Bartos et al. 
(2015) suggested that on-tlie-fly Ha mapping may 
be used to generate catalogs of galaxies for kilo- 
nova searches triggered by GW events, and found 
that mapping a ~ 10 deg 2 region to identify these 
galaxies within 200 Mpc can be done cost effec¬ 
tively and quickly, within a day, using a 1-2 meter 
class telescope. 

Finally, another approach for a targeted sur¬ 
vey is to focus on the galaxies in the 50% (as op- 


5 



posed to 90%) confidence region of the localized 
sky area associated with a GW event. Such an 
approach can decrease the covered sky area by a 
factor of ^5 (Singer et al. 2014). In this case, if 
the targeted survey follows up three GW triggers 
instead of one, then a similar probability of success 
(90%) is achieved by covering 40% fewer galaxies 
than the survey of a single 90% confidence region. 
Such an approach, however, does not require 40% 
less time since the slew to the first galaxy follow¬ 
ing a GW trigger would occur three times. Tar¬ 
geted surveys of all nearest (<200Mpc) galaxies 
and only nearby Ha galaxies within 50% confi¬ 
dence regions of the localized areas of three op¬ 
timal GW sources could be done with NIRCam 
within a total of 8.7 hours and 4.4 hours, respec¬ 
tively. Not only would such an approach require 
less time, but it would also result in fewer galaxies 
needing followup observations to confirm a kilo- 
nova identification. 

JWST is not an efficient facility for a survey 
involving quick exposures of many fields. In these 
examples of targeted surveys, the overhead times 
are 50-70 times greater than the exposure times. 
Despite this inefficiency, by judiciously choosing 
the galaxies most likely to host kilonovae associ¬ 
ated with GW events, the NIRCam observing time 
can be minimized to enable high-impact science 
not possible from other facilities. 

5. Identification of kilonovae 

Kilonovae and supernovae are among the intrin¬ 
sically brightest extragalactic compact sources, 
more than two orders of magnitude brighter than 
other, more typical sources at optical wavelengths 
(e.g., Kasliwal 2013; Rau et al. 2009). In the near- 
infrared, kilonovae are even brighter (Barnes & 
Kasen 2013). A single-epoch NIRCam observa¬ 
tion is therefore sufficient to identify any galaxy 
associated with a kilonova candidate, particularly 
a candidate identified within a week of and coin¬ 
cident with a GW trigger event. 

Once a kilonova candidate is identified, one 
needs to distinguish it from the foreground (e.g., 
Milky Way stars and asteroids), background (e.g., 
distant galaxies), and unresolved (e.g., HII re¬ 
gions) sources. For example, multiple background 
galaxies that are sufficiently bright may overlap 
a nearby (~ 200 Mpc) galaxy, resulting in a con¬ 


tinuous quasi-point source with brightness com¬ 
parable to kilonovae. Distinguishing between a 
bona-fide kilonova and continuous sources requires 
a template be obtained, either before or after the 
occurrence of the kilonova, or multiple NIRCam 
exposures obtained days apart. For the former, 
there is currently no suitable all-sky survey for 
this purpose, but it will be possible to carry out a 
survey selectively for the relevant galaxies, using 
other facilities to complement NIRCam observa¬ 
tions. For the latter, NIRCam followup observa¬ 
tions could be obtained on the galaxies. This ap¬ 
proach may not be optimal since it requires double 
JWST time, but the advantage is that such obser¬ 
vations would provide the template to distinguish 
kilonova from continuous sources and enable char¬ 
acterization of its dimming. 

Some transient sources (e.g., foreground aster¬ 
oids and dwarf novae, background supernovae) 
could be mistaken as kilonovae, resulting in false 
positive transient events. To address this issue, 
we first recall the results of Kulkarni & Kasliwal 
(2009), which finds that the background rate of su¬ 
pernovae, aligned by chance with a nearby galaxy, 
for Advanced LIGO-Virgo in a single snapshot 
within 12 deg 2 with r-band luminosity of r < 24 
should be < 0.1. That study also finds the fore¬ 
ground rate of flares to be even less, < 0.01. These 
sources are expected to be the dominant source 
of false positive events. While these false positive 
event rates could, in principle, be different in NIR- 
Cam’s infrared wavelength range, they are unlikely 
to be significantly greater. We therefore expect 
that any survey focused on galaxies, and especially 
the targeted surveys discussed in §4, will render 
the number of false positive transient events prac¬ 
tically negligible. 

6. Probing the kilonova emission model 

The detection of kilonovae from a compact bi¬ 
nary merger detected via GWs will help answer 
a number of important questions. Many of these 
questions may be addressed with the detailed ob¬ 
servation of one kilonova, making even the lim¬ 
ited observation time required for one detection 
valuable. If one determines the peak flux and the 
timescale for the kilonova to reach maximum light, 
one can deduce the quantity of r-process ejecta and 
its mean velocity (e.g., Metzger et al. 2010). For a 
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statistical sample of such events, this information 
will address whether neutron star mergers are the 
dominant source for producing r-process elements 
in our Galaxy. When coupled with GW measure¬ 
ments of the parameters of the merging binary, 
the mass and velocity of the ejecta will also help 
constrain the equation of state of dense nuclear 
matter (e.g., Bauswein et al. 2013). 

Determining the peak flux of a kilonova may 
nevertheless not be possible unless the light curve 
is observed multiple times during the emission pe¬ 
riod. In principle, such rapid follow-up observa¬ 
tions are possible. The initial survey, executed 
as a target of opportunity, could occur within 
2 days (operations requirement MO-210; Jordan 
2014) of the GW event. With calibrated, pro¬ 
cessed NIRCarn images then available within 5 
days after downlink (e.g., operations requirement 
MO-41; Jordan 2014), follow-up NIRCarn obser¬ 
vations of a candidate or confirmed kilonova could 
then be executed as a target of opportunity again 
within 2 days. Thus, if an investigator develops 
data processing pipelines such that an appropri¬ 
ate gravitational wave event can be selected and 
NIRCarn images be searched in relatively short pe¬ 
riods of time, two epochs of NIRCarn observations 
separated by a week may be obtained for a kilo¬ 
nova within 9-10 days of the binary merger. This 
time frame is likely sufficiently short to enable ad¬ 
ditional observations of the kilonova, providing in¬ 
formation on its temporal evolution. 

Triggered spectroscopic follow-up, either with 
JWST itself or a large ground-based near-infrared 
telescope (e.g., the Giant Magellan Telescope), 
could confirm the merger origin of the event by 
detecting the absorption lines of exotic r-process 
elements (e.g., Kasen et al. 2015b). The strength 
of individual lines, once identified, could in princi¬ 
ple be used to determine the relative abundances 
of individual nuclei. 

7. Conclusion 

We showed that JWST/NIRCarn can easily de¬ 
tect kilonovae out to distances relevant to GW 
observations of compact binary mergers. To ef¬ 
ficiently survey the sky for kilonovae following a 
GW detection, NIRCarn observations will need to 
be directed toward galaxies within the GW dis¬ 
tance range. For a maximum source distance of 


200 Mpc, which will be the typical distance of de¬ 
tected binaries with Advanced LIGO-Virgo at de¬ 
sign sensitivity, the required NIRCarn observation 
time for a fiducial 10 deg 2 GW-localized sky area 
would be about 13 hours overwhelmingly domi¬ 
nated by slewing and guide star acquisition. This 
time can be significantly decreased, perhaps to 
about 4 hours, by focusing on (i) the galaxies with 
the highest star formation rates and (ii) especially 
those within the most probable GW sky area, en¬ 
abling a survey of multiple GW detections in less 
time without compromising detection of a kilonova 
counterpart. 

We find that the identification of kilonovae is 
unlikely to be limited by foreground or background 
transient events. With such detection capabil¬ 
ity, JWST/NIRCam will be able to (i) regularly 
survey GW event candidates for kilonova counter¬ 
parts, therefore establishing a statistically signifi¬ 
cant sample size of kilonova emission parameters, 
and (ii) will be able to rapidly find kilonovae, en¬ 
abling detailed study of their temporal and spec¬ 
tral evolution. 
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